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1 Introduction 



Vector meson electroproduction at large photon virtuality, Q^, has always attracted a lot 
of theoretical interest. Its diffractive nature as well as the interesting correlation between 
the and the energy dependence are challenging issues. At first traditional concepts like 
vector meson dominance (see e.g. [Tj) or the Regge model with its prominent Pomeron 
exchange (see e.g. j2]) have been exploited to analyze the electroproduction data. In 1987 
Donnachie and Landshoff 3J viewed the Pomeron as the exchange of two gluons between the 
proton and the quark- ant iquark pair created by the virtual photon and which subsequently 
form the outgoing meson. Brodsky et al. jl] treated the two-gluon exchange contribution 
to electroproduction at large and small Bjorken-x, xbj, in the framework of QCD 
factorization. They showed that in their approach, known as the ln(l/a;Bj) approximation, 
the emission and reabsorption of the gluons by the proton can be related to the usual 
gluon distribution. Many variants of the leading ln(l/xBj) approximation can be found 
in the literature which differ mainly by the treatment of the subprocess 7*(? — > Vg, see 
Ref. 13 El 13 IHI to name a few. These approaches describe many features of vector meson 
electroproduction quite well. 

In 1996 vector meson electroproduction has been taken up by theory again. Exploiting 
the new concept of generalized parton distributions (GPD) jHI El it has been shown |H1 
ITT| that, at large Q^, the process factorizes into a hard parton- level subprocess - meson 
electroproduction off partons - and soft proton matrix elements which represent generalized 
parton distributions. The process is depicted in Fig. ^ where also the momenta of the 
involved particles are specified. It has also been shown in Refs. E] that the process is 
dominated by transitions from longitudinally polarized photons to longitudinally polarized 
vector mesons (L — > L) at large Q^] the amplitudes for other transitions are suppressed by 
inverse powers of Q. The production of vector mesons at small xbj ( ^ 10^^) is controlled 
by gluonic GPDs where quasi on-shell gluons are emitted and reabsorbed by the proton. 
These GPDs which represent the soft physics embodied in the proton matrix elements, are 
unknown as yet and have to be modelled. 

Detailed experimental information on electroproduction of light vector mesons in the 
region of small xbj is available from HERA. Cross sections and spin density matrix elements 
have been measured by HI [T21 and ZEUS [121 Cll. Despite the sound theoretical basis 
of the handbag approach not much has been done as yet in analyzing these data within 
this framework. There is only the explorative study of the longitudinal cross section for p 
production performed by Mankiewicz et al. ^^1- The normalization of the cross section was 
however not understood in this work. Martin et al. on the other hand, started from 
the ln(l/xBj) approximation and estimated effects due to the replacement of the gluon 
distribution by the corresponding GPD. Here, in this work we attempt a complete and 
systematic analysis of the available electroproduction data at small x-Qy In order to analyze 
the spin density matrix elements of the vector mesons we also calculate the amplitudes for 
transitions from transversely polarized photons to transversely and longitudinally polarized 
vector mesons (T ^ T and T L). We allow for quark transverse momentum and take 
into account Sudakov suppressions. As it will turn out this approach leads to the correct 
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Figure 1: The handbag-type diagram for meson electroproduction off protons. The large 
blob represents a GPD while the small one stands for meson electroproduction off partons. 
The momenta of the involved particles are specified. 

normalization of the cross sections at finite but large . Infrared singularities which occur 
for the T ^ T transition amplitude in coUinear approximation ^Hj, are also regularized in 
our approach although in an admittedly model-dependent way. 

The plan of the paper is the following: A kinematical prelude and the handbag ampli- 
tude are presented in Sect. 2. The amplitudes for the subprocess — >• Vg are discussed 
in Sect. 3 to leading order of perturbative QCD and including transverse momenta of the 
quarks and antiquarks making up the meson. The impact parameter representations of 
the full handbag amplitudes for electroproduction of vector mesons are presented in Sect. 
4. The following section, Sect. 5, is devoted to the construction of the GPDs. Numerical 
results, obtained from the handbag approach, for the cross sections of vector-meson elec- 
troproduction and for the vector meson's spin density matrix elements are compared to 
recent experimental results in the small xbj region in Sects. 6 and 7, respectively. In the 
next section. Sect. 8, we discuss the helicity correlation Ai^i^ and the role of the GPD H 
and summarize in Sect. 9. 



2 The handbag factorization 

We will work in a photon-proton center of mass (cm.) frame, see Fig. Q in a kinematical 
situation where 

W = {v^ qf , (1) 
and the virtuality of the incoming photon, = —Q^, are large while Bjorken's variable, 

XBj = QV{2p-q), (2) 

is small (xbj ^ 10^^). We also assume the square of the momentum transfer, A = p' — p, 
to be much smaller than Q^. 
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In light-cone components, defined by = (a°±a^)/-\/2 and a-b = a^h +a 6+ — a^bj 
the momenta of the protons and the photon read 
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(3) 



where t] equals — xbj up to corrections of order w? jQ"^ and A^/Q^. Here, m denotes the 
mass of the proton. The average proton momentum is defined by 



P 



and the skewness parameter ^ by 



2 (P + P ) ' 



{p - p'y 



(4) 



(5) 



In the photon-proton cm. frame and for small xbj, the skewness parameter is related to 
Bjorken-x by 



XBj/2. 



For Mandelstam t, given by 



t 



2-a;Bj 

AO^rri^ + A^ 



1-e 

a minimal value is implied by the positivity of 

e 



(6) 



(7) 



4m^ 



Since we are interested in the region of small Bjorken-x and, hence, small skewness we will 
use train — in the following. We also will neglect the proton and meson (my) masses in 
the kinematics. 

Let us now consider the dynamics of vector meson electroproduction in the kinematical 
regime specified above. The dominant contribution in this kinematical region comes from 
the emission and reabsorption of collinear gluons from the protons accompanied by '-y*g —>■ 
Vg scattering [H]. The neglect of an analogous quark contribution is justified by the 
fact that, at small xbj, partons with small momentum fractions dominantly participate in 
hard meson electroproduction. Since, at small — t, the GPDs are expected to reflect the 
magnitudes of the usual parton distributions the gluon contribution should outweigh the 
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quark one. This is in particular the case for electroproduction of mesons where only the 
small strange quark content of the proton is probed. Even for the production of p mesons 
the gluonic contribution seems to be still sizeable for xbj as large as 0.1 as is indicated by 
the ratio of and p electroproduction cross sections [T7] . 

The momenta of the gluons which, approximately, are coUinearly emitted or absorbed 
by the protons, are parameterized as 



kn 



k' 



8(x + Or 



-,-A^/2 



8(x - ^)p+ 



A^/2 



In general the partons may have small virtualities of the order of 
introduced an average parton momentum 



(9) 



As usual we have 



k 



-{kg + k' 



(10) 



and an average momentum fraction 



X 



k+/p-^ 



(11) 



In order to facilitate comparison with other work we also provide the relations between 
the variables x and ^ and the usual Mandelstam variables for the hard subprocess. They 
read (t ^ 0) 



u 



{q + kg 
{q-k 



I \2 



2i 



(12) 



and are valid at large and small ^. 

Radyushkin has calculated the asymptotically leading handbag contribution to me- 
son electroproduction at small xbj |9;. As he showed this contribution involves L — > L 
transitions. Leaving aside for the time being a potential breakdown of factorization, 
Radyushkin's result can straightforwardly be generalized to other transitions fHl ^] • The 
crucial point in the derivation of the handbag amplitude is the use of light-cone gauge for 
the gluon field, n ■ A" = 0, where 



n 



[0,1, Oj 



(13) 



and a is a color label. This gauge allows to express the gluon field by an integral over the 
gluon field strength tensor Gl^i jHl HH] (the limit e — is to be understood) 



n 



dae-''' Gt,Az + an) 



(14) 



5 



With the help of this expression one can replace the products of fields appearing in the 
perturbatively calculated amplitude for ~^ by 

raa' „ 
^^c ^ A,A'=±1 

X n-' GlJa'n) G^^iz + an) e*''{kg, A) e^' {k'^, A') , (15) 

where we have also made a helicity projection for the gluons. The use of the approximation 
Q for the gluon momenta forces the relative distance of the fields on the light cone 
z z = [0,2;_,0_l]. The vectors e{kg,X) and e(fc^, A') specify the polarization of the (on- 
shell) gluons, the corresponding momenta, kg and k'^, are defined in Eq. ©. The first set 
of polarization vectors in (jl5|) is to be used to contract the hard scattering kernel leading 
to gauge invariant parton-level helicity amplitudes Ti-^iy^fj^x for Vg (/x and /x' denote 

the helicities of 7* and V, respectively). The contraction of the field strength tensors with 
the second set of polarization vectors leads to 



n-n-' GU^^'n) G,.^.{z + an) e*'^(A;„ A) e^'(fc;. A') = 

In^n-' G,^{a'n)G,,A'z + an) [ (-^^f + A^ef ) 5aa' - ^a-vI , (16) 
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where 



„ii _ „22 _ _ 12 _ 21 _ ,11 _ ,22 
y± — y± — ^± — ^± — — 
.12 J.21 



t]^ = ei = i\ (17) 

while all other components of these tensors are zero. That only the transverse components 
in the contraction remain is a consequence of the chosen light-cone gauge and of the fact 
that the polarization vectors have zero plus components in the cm. frame we are working. 

Proton matrix elements of the gluon helicity non-flip contributions and ie^^ in 
(fTBj) define the unpolarized, H^{x,^,t) and E^{x,$,,t), and the polarized, H^{x,^,t) and 
E^{x,^,t), gluon GPDs, respectively [Hllini- The proton matrix elements of these gluon 
field operators are related to the GPDs by 

{p'u'\Y^A^^iO)A-'^'iz)\pu) = 1 5: e^{kg,\)e*'''{k'g,\') (18) 

a, a' A=±l 

Jo {x + C, — ie){x — + ie) 
y 2p ■ n Am p ■ n 

2p ■ n Am p ■ n J 

Working out the spinor products one sees that for proton helicity non-flip the linear com- 
binations [20] 

H^{x,i.t)-^L_E3{x,i,t) (19) 
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and 



(20) 



occur. Since we are interested in small ^ the and terms can safely be neglected in 
the expressions and ()2()j) . For proton helicity flip, on the other hand, and do 
not contribute but only 

~^ ^ E^2u^E'^), (21) 



— K 



2m 1-C 



where k is a phase factor reading 



(22) 



for proton momenta of the form (jS)). 

The gluon helicity flip contribution in (fTBj) which defines four more GPDs |2I], will be 
neglected in the following since it is strongly suppressed at small —t. The mismatch between 
the proton and gluon helicities in the proton matrix elements has to be compensated by 
orbital angular momentum. For each unit of it a factor y/^/m is picked up fTI\ 1^ . 
Further suppression comes from the subprocess amplitudes which behave as 



njV 



(23) 



at small —t and from the fact that the amplitude Hq-a^a vanishes for /i = ±1 [H 

Combining all this, we finally obtain the helicity amplitudes for electroproduction of 



vector mesons 



3. 
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(24) 



for proton helicity non-flip (explicit helicities are labelled by their signs) and for helicity 
flip 



+ 



-- LyK 

2 2m 



dx 



X + ^){x - ^ + ie) 
E\x,^,t) 

^E%x,U)} ■ 



(25) 



•^We note in passing that our helicities are hght-cone hehcities which naturally occur in the handbag 
approach. The difference to the usual cm. frame helicities is of order m^/^/W'^ 21 and can be ignored 
in the kinematical region of interest in this work. 
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Figure 2: Lowest order Feynman graphs for the subprocess 7*(7 Vg. 

are functions of Q^, x, ^ and t. The flavor weight factors, 

-ea) = 1/V2- C^ = es = -1/3. (26) 

where denotes the quark charge in units of the positron charge e. The remaining hehcity 
amphtudes are obtained with the help of parity invariance 

M-,,-,,,-,-. = (-1)^'— M^,.,,^.. (27) 

An analogous relation holds for the subprocess amplitudes. 

3 The partonic subprocess 7*^ Vg 

The parton-level amplitudes for the subprocess 7*(7 Vg are calculated from the Feyn- 
man graphs shown in Fig. El the outgoing qq pair is to be combined into the vector meson 
regarding its quantum numbers. This is conveniently done by means of a covariant spin 
wavefunction. As is well-known from analyses of hadron form factors at large momentum 
transfer, leading-twist perturbative calculations are instable in the end-point regions since 
the contributions from large transverse separations, b, of quark and antiquark forming 
the meson are not sufficiently suppressed. In order to eliminate that defect Li and Ster- 
man retained the quark transverse degrees of freedoms and took into account Sudakov 
suppressions. Including, in addition, meson wavefunctions with their intrinsic transverse 
momentum dependence instead of distribution amplitudes |21], the perturbative contribu- 
tions to form factors can reliably and self-consistently be calculated, the end-point regions 
are strongly damped. 

Since the subprocess 7*5' — » Vg bears resemblance to the meson form factors it is 
tempting to apply this so-called modified perturbative approach also here in order to 



The subprocess amplitudes, 7i , 
Cy, read for p and mesons 



1 

3 
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suppress the contributions from the soft end-point regions and, simultaneously, to regularize 
this way infrared divergencies that may occur in the T L and T ^ T amplitudes 
[T6] . The modified perturbative approach applied to the subprocess, is, to some extent, 
similar to the mechanism proposed in Ref. 6j for the suppression of the leading-twist gluon 
contribution to hard meson electroproduction. It is however to be stressed that in Ref. [H] 
the leading ln(l/xBj) approximation of Ref. |1] has been utilized. 

Let us now turn to the description of the soft qq V transtion matrix element. We 
start from a frame where the hadron rapidly moves along the 3-direction (g' = [0, q'~ , 0±] 
with q'~ ~ Q/(2v^)). This frame is termed the hadron-out one. The momenta of the quark 
and the antiquark which form the valence Fock state of the meson, are parameterized as 

q'i = Tq"^ + K, qti = fq'^ + k^. (28) 

where 

fci = [A:+,0,kiJ, k2 = [kt,Q,k2±]- (29) 

The variables r and f are the usual fractions of the light-cone minus- component of the me- 
son's momentum the constituents carry. Momentum conservation provides the constraints 

f = 1 - r , k2± = -ki^ = -k^ . (30) 

It can be shown (23] that the variables r, f and k^ are invariant under all kinematical 
Poincare transformations, i.e. under boosts along and rotations around the 3-direction as 
well as under transverse boosts. Moreover - and this is an important point - the light- 
cone wavefunction associated with the valence Fock state, \l/y = \l/i/(r, kj_), is independent 
of the hadron's momentum and is invariant under these kinematical transformations too. 
The light-cone wavefunction may differ for longitudinally and transversally polarized vector 
mesons j2S|- 

As is customary in the parton approach we neglect the binding energy. That possibly 
crude approximation can be achieved by putting the individual fc^-components to zero. In 
fact starting from a parameterisation of the various momenta in the meson's rest frame and 
boosting to the hadron-out frame, one sees that the fc^-components are of order my/q'~. 
The plus-component of the difference of the momenta 

K =^ik,- k,) , (31) 

is zero with this choice of k^ components and, hence, 

K = [0, 0, kj , (32) 

and q' ■ K = 0. The quarks are treated as massless in the hadron-out frame; they are not 
strictly on- shell. 

It is convenient to couple the spinors representing quark and antiquark in a covariant 
spin wave function for the vector meson. The Dirac indices of it (omitted for convenience) 
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are to be contracted with the corresponding ones of the hard scattering kernel (see below) . 
For the construction of the spin wave function we adapt the method presented in Ref. j2ZI 
(see also .28j) straightforwardly to vector mesons. The product of spinors v{q2)u{qi) is 
boosted to the hadron's rest frame, coupled there into the quantum numbers of the vector 
meson and boosted back to hadron-out frame. Separating terms with and without K and 
neglecting terms oc , one arrives at 

Tv = Tvo + AFy^iT" , (33) 

where ^ ^ 

r^o = {i' + mv)liv, Ar„ = — {ryo,7a} • (34) 

The polarization state of the meson is described by the vector e. The soft physics parameter 
My is of order mv'-, its model dependence results from the specific treatment of the quarks 
in the meson's rest frame. In the following we will use My = my for simplicity but we will 
comment on other choices of it. 

Since the anticommutator {Fyo; i'} is zero the 4-vector K is only determined up to a 
multiple of the meson momentum. This property can be used to identify K^, given in Eq. 
(I32|) . with the quark-antiquark relative momentum 

^ ^-{q^ - q,) , (35) 

where the parton momenta, qt, are defined in Eq. (j28|l . This choice, although not forced, 
is very convenient. Its main advantage is that, K'^ now represents one unit of orbital 
angular momentum in a covariant manner j^Zj. As discussed in this article, the relative 
momentum ()35p is a 4-transverse vector which are defined by Kj^ = K'^ — q' ■ K/myq'^. 
In the hadron-out frame and up to corrections of order my/q'~ ., K± = K. In the meson's 
rest frame on the other hand, clearly K_[^ = (0, k), and one has an object transforming as 
a 3- vector under the three-dimensional rotation group 0(3). 

One of the basis ingredients of the hard scattering picture is the collinear approximation 
which says that all constituents move along the same direction as their parent hadron up 
to a scale of the order of the Fermi motion (k^) which typically amounts to a few 100 MeV. 
The (nearly) collinear kinematics justifies an expansion of the spin wavefunction upon a 
power series in or, in order to retain a covariant formulation, in K^. Up to terms linear 
in K^^ this expansion is given above for vector mesons. 

The transformation from the hadron-out frame to our cm. frame where the meson 
momentum has a transverse component — A_|_/2, is executed by a transverse boost (cf. 
e-g- ES) that leaves the minus component of any momentum vector a unchanged, and 
which involves a parameter d~ and a transverse vector d_|_, is defined as 

r + - 1 r + a_L ■ d_L a~ d| _ j i /q«\ 

[a^,a ,a±\ — ^ [a^ — h ^ , a , a_L - ^ d^J . (36) 

The transverse boost is one of the kinematical Poincare transforms that leaves the hadron 
wavefunction invariant. Taking for the parameters d~ = q'~ and d_L = A_l/2, we readily 
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find from ()28p tlie expressions for the quark momenta and the relative momentum in the 
cm. frame. 

Provided quark transverse momenta are taken into account, the general structure of 
the 7*5^ Vg amphtude is 

= / ^vir, kl) Tr [r^ U . (37) 

The hard scattering kernel, Th = Th{t, x, Q^, K, t), can be written as follows 

Th = To(r, x, Q\ k\, t) + AT^(r, x, Q\ k\, t)K^ + . . . . (38) 

Terms oc K^K^ (and higher) in the nominator are neglected while, in the spirit of the 
modified perturbative approach [23], the k\ terms in the denominators are kept. Inserting 
Eq. as well as the spin wavefunction (j3^ into Eq. (jHTjl . one obtains 

vl>^(r, kl) Tr [Tvo T, + ATy^ T, 
+ro AT^ + AEy, A iT" + . . .] . (39) 

Obviously, the terms oc integrate to zero while the K^K'^ term survives the k±_- 
integration . Hence, 

= / dT^^^>y{T,kl)l:T^^Ty,T,-hlgf ATvc.^Tp + ..^ , (40) 

where g±_ is the transverse metric tensor defined in Eq. (|T7j) . In order to simplify matters 
we only take into account the first non-zero term in this expansion for each amplitude, 
i.e. we neglect any correction of order my or k^ to its leading term ^. As we said above 
we however retain the kl terms in the denominators of the propagators. Morever, any t 
dependence of the subprocess amplitudes is ignored except the factors of required by 
angular momentum conservation. This is justified in the small t region we are interested 
in. 

For longitudinally polarized vector mesons the first term in Eq. (j40p contributes, the 
other term represents a kl correction to it which we, according to our strategy, neglect as 
well as all other terms indicated by the ellipses. For transversely polarized mesons, on the 
other hand, the first term in Eq. ()40|) disappears since the number of 7 matrices in the 
trace is odd ^. The second term in Eq. (jlOI), oc k\, contributes in this case; it scales as 
oc k\/ {iriyQ), see Eq. ()34|1 . Combining this property with the behavior of the subprocess 
amplitudes near the forward direction (^3)) and utilizing the fact that ik"])^^"^ /my is of 
order 1, the various photon- meson transitions respect the following hierarchy 



"^Note that the hard scattering kernel Th does not depend on the vector meson mass; it occurs through 
the spin wave function. 

^We remind the reader that for longitudinaUy polarized vector mesons e(0) = q' /niy up to corrections 
of order my /q'~ ■ 
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L ^ L : T^oxfi A 0^ 1 5 



T ^ L : ^^A,+A 



T-^T: n\^,^ 



+A,+A ^ /-) ) 

L ^ T : ^+A,oA 



Q ' 



T^-T: „v.^,,,.^<^±2_. (41) 

This hierarchy propagates to the proton non-flip amphtudes for the full process and justifies 
the neglect of L — > T and T —>■ —T transitions in the analysis. The amplitudes for proton 
helicity flip have an extra factor i/^/m, see (EHjl . Our interest in this work is focussed on 
unpolarized protons. In the corresponding cross sections there is no interference between 
flip and non-flip amplitudes. Hence, proton flip is suppressed by a factor of t and since 
there is no theoretical or phenomenological indication that \E^\ is much larger than , 
neglected by us. Information on the proton flip amplitudes may be extracted from data 
on meson electroproduction with polarized protons. As a last simplification we neglect 
contributions from in the evaluation of the amplitudes. Since in the forward limit 

— > 0, and reduce to xg{x) and xAg{x), respectively, it is plausible to expect 
that the relative size of Ag and g is reflected in that of and at small C, and —t. 
Since |A(7(x)| is much smaller than g{x) the contribution from can safely be neglected. 
The model GPDs we are going to construct in Sect. Eldo indeed respect this assertion. As 
a consequence of parity invariance, see ()24|) and (j^7|) . there is anyway no contribution from 
the GPD to the most important amplitude, L ^ L. Care is required for observables 
for which the contribution from partially if not totally cancels. An example of such 
an observable is the correlation of the electron and proton helicities. We will comment on 
this observable in Sect. |H1 

The hard scattering amplitudes for the three helicity configurations we keep in our 
analysis are to be calculated from the Feynman graphs shown in Fig. |21 The results for 
the relevant sums and differences of positive and negative gluon helicities can be cast into 
the following form 

K... ± K-.- = i'* / (^''e)fp,D, (42) 

where the product of propagator denominators reads 

>i -fix-O QV(2 - ^e) {kl + f (x + QVi2 
X (kl + r{x + 0QV{20) (kl-Tix-0QV{20-^e) . (43) 



12 



Here, Nc denotes the number of colors. The functions /^^ read 
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kl 




Tf 
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mv 
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kl 
my 
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/oo 


— fo+ - 


-- 0. 







2^3 



(44) 

Following Ref. we retain fc^ terms in the denominators of the propagators ()43p . These 
terms play an important role since they compete with terms oc T{f)Q^ which become small 
in the end-point regions where either r or f tends to zero. They lead to the suppression 
of contributions with large quark- ant iquark separations as we mentioned above. 

In collinear approximation and utilizing distribution amplitudes up to twist-3 accuracy 
the subprocess amplitudes for T — T transitions are infrared divergent, signaling the break 
down of factorization ^H]- Neglecting transverse momenta in Eq. ()43|) . one finds 

+ "'"^'j/^y / *i *i *vT(r, k'i) . (45) 

Assuming for instance a Gaussian wavefunction "^vt ~ exp [(Ivt^±/{'^^)]^ ansatz that 
has been shown to work successfully in many cases (see for instance [23]) and will be used 
by us in the numerical analysis of meson electroproduction, we find that in fact the r 
integral is regular. As pointed out in [WJ, the x integral in Eq. (j24j) may not exist due 
to the double pole (x — ^ + ie)"^ occurring. Whether or not this happens depends on 
properties of the GPDs. In Sect. El we will take up this problem again. 

One may also consider a transverse momentum dependence of the GPDs. That issue 
has been investigated in Ref. POj for meson electroproduction at intermediate values of 
XBj. In this kinematical region the emission and reabsorption of quarks from the proton 
dominates. We however think that the k± dependence of the GPDs is of minor importance. 
In contrast to the meson where the hard process enforces the dominance of the compact 
valence Fock state of the meson, all proton Fock states contribute to the GPDs at small 
—t \20\ 1^ . If the gluons are distributed in the proton like the quarks, an assumption 
that is supported by the slope of the diffraction peak in elastic proton-proton scattering, 
the k±_ dependence of the GPD should roughly reflect the charge radius of the proton 
((^1)^'^^ — 200 MeV). Consequently, we expect to be only mildly dependent on the 
transverse momentum, a potential effect we neglect. 
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4 The impact parameter space 



Transverse momenta in the subprocess amplitudes, see Sect. El imply finite quark-antiquark 
separations in the configuration space which are accompanied by gluon radiation. On the 
grounds of previous work by Collins and Soper Sterman and collaborators |2S] cal- 
culated this radiation to next-to-leading-log approximation using resummation techniques 
and having recours to the renormalization group. The result is a Sudakov factor which 
suppresses large quark-antiquark separations and which we also have to take into account 
in our analysis in order to have consistency with the retention of the transverse degrees of 
freedom. Since the Sudakov factor is given in the transverse separation or impact param- 
eter space - only in this space the gluonic radiative corrections exponentiate - we have to 
work in this space. 

The two-dimensional Fourier transformation between the canonical conjugated b and 
k± spaces is defined by 

/(b) = J d'k^ exp Hk^ • b] /(k^) . (46) 

For the meson wavefunctions we adopt the same Gaussian parameterization as is used for 
the pion PH E21 



TT 



(47) 



{i = L, T) which strictly speaking represent full wavefunctions with their perturbative tails 
removed. Transverse momentum integration of these wavefunctions lead to the associated 
distribution amplitudes which represent the soft hadronic matrix elements entering calcu- 
lations within the coUinear factorization approach. Actually, the wavefunction ()47|1 lead 
to the so-called asymptotic meson distribution amplitude 

$As = Qrf . (48) 

For the decay constants, fvi of longitudinally polarized vector mesons we take the values 



fpL = 0.216 GeV, Ul = 0.237 GeV . (49) 

The decay constants for transversely polarized vector mesons are almost unknown. The 
only available information comes from QCD sum rules. In Ref. j^H] fpT has been estimated 
to (160 ± 10) MeV. We actually fit these decay constants to experiment. Identifying for 
instance the soft parameter My in the spin wave function with the meson mass, we obtain 

fpT = 0.250 GeV, /<^t = 0.275 GeV. (50) 

Choosing My to be smaller than the meson mass results in smaller values of the decay 
constants fvT- The transverse size parameters ayi are fixed by the requirement of equal 
probabilities for the vector meson and pion valence Fock states, namely 0.25. This leads 
to 

apL = 0.52 GeV-\ a^L = 0.45 GeV'^ (51) 
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The transverse size parameters for transversely polarized vector mesons are adjusted to 
experiment. The numerical results we are going to present below are obtained with 



apT = 0.65 GeV"^ , a^r = 0.60 GeV^^ . 



(52) 



With the parameter values quoted in and ()5H1 the r.m.s. transverse momenta, evalu- 
ated from (jUj), amount to 0.61 GeV and 0.67 GeV for the longitudinally polarized p and 
mesons, respectively. These values are much larger than the one for the proton r.m.s. 
transverse momentum. 

The Fourier transform of the meson wavefunction (jTfjl reads 



-rr 



(53) 



The product of the propagator denominators D (|l3|) can be decomposed into single-pole 
terms which are either of the form 



1 



or 



1 



ki - d2{x - OQ^ - le ■ 
where di > 0. The Fourier transforms of these pole terms can readily be obtained 

1 



(54) 
(55) 



Ti 



27T 
1 



T2 = —Ko{^d2{^-x)bQ) e{^~x) 



+ '-H^''>{^d2{x-0bQ] 9{x-0 



(56) 



where Kq and Hq^^ are the zeroth order modified Bessel function of second kind and Hankel 
function, respectively. 

Putting all this together and including the Sudakov factor, exp[— S'(r, b, Q)], the gluonic 
contributions to the helicity amplitudes for vector-meson electroproduction read 



Cv JdxdT f+^H<^{x,U) 



V2K 

X / d%^v,.'ir,b'^)D{T,Q,b)as{pR) exp [-5(r, 6, Q)] 
Cv / dxdr f-, H<^{x,^,t) 



d% ^v^.'{r, b^) D{r, Q, b) exp [-S{t, b, Q)] 



(57) 
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which is the b-space version of the amphtude ()24|) . The functions D and fp^ are given 
in PHj) and (jH}. Since the Fourier transformed wavefunctions, the product of propagator 
denominators as well as the Sudakov factor only depend on b the angle integration in the 
last integral is trivial and a three-dimensional integral [dxdrb dh) remains to be evaluated 
numerically. Parity invariance ^TJ\ lead to the following relations among the amplitudes ^ 

-^1+,++ = , Mt,+^ = Mt,-+ , (58) 

The Sudakov exponent S in (j^Zj) is given by ^2S1 

5(r, 6, Q) = s{t, b, Q) + 3(f, 6, Q) - ^ In M^^QCd) ^ ^^^^ 

where a Sudakov function s occurs for each quark line entering the meson and the abbre- 
viation 

6 = -ln(6AQCD), (60) 

is used. The last term in ()59|) arises from the application of the renormalization group 
equation {(3q = 11 — |n/) where nj is the number of active flavors taken to be 3. A value of 
220 MeV for Aqcd is used here and in the evaluation of as from the one- loop expression. 
The renormalization scale fiR is taken to be the largest mass scale appearing in the hard 
scattering amplitude, i.e. fiR = max (rQ, fQ, 1/6). For small b there is no suppression 
from the Sudakov factor; as b increases the Sudakov factor decreases, reaching zero at 
b = 1/Aqcd- For even larger b the Sudakov is set to zero The Sudakov function s reads 

s(r, b,Q) = ^ (^q In (fj - q + b^ + NIL - terms , (61) 

where 

q = In (rQ/( V2Aqcd)) • (62) 

Actually we do not use the explicit form of the next-to-leading-log corrections quoted in 
but those given in Ref. |3S1- The latter ones contain some minor corrections which 
are hardly relevant numerically. Due to the properties of the Sudakov factor any contri- 
bution to the amplitudes is damped asymptotically, i.e. for lia{Q^ / Aqq^^) oo, except 
those from configurations with small quark-antiquark separations, b plays the role of an 
infrared cut-off; it sets up the interface between non-perturbative soft gluon contributions 
- still contained in the hadronic wavefunction - and perturbative soft gluon contributions 
accounted for by the Sudakov factor. 

^The same relations as for the terms also hold for the t-channel exchange of a particle with natural 
parity, P = (—1)'^. The relations for the terms in (|57|l corresponds to those obtained for an unnatural 
parity exchange |34j . 

^The definition of the Sudakov factor is completed by the following rules [221: exp [—5*] = 1 if exp [—5*] > 
1, exp [-5"] = if 6 > 1/Aqcd and s(/3, 6, Q) = if 6 < V2/PQ. 
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5 Modelling the GPDs 



In order to calculate the electroproduction amplitudes (jSlJ we still need the gluon GPDs. 
A model for them can be constructed with the help of double distributions [HE] which 
guarantee polynomiality of the GPDs. The gluonic double distribution f{P,a,t ~ 0) is 
customarily parameterized as 

g t ~ 0) = „(m ^P" + ^> I'l-W-^T ,63) 



where g{x) is the usual gluon distribution. Its definition is extended to negative (3 by 

9{-P) = -9{f3) . (64) 

A popular choice of is 1 for quarks and 2 for gluons. This is motivated by the interpreta- 
tion of the a dependence like a meson distribution amplitude for hard exclusive processes. 
The cases n = 1 and 2 correspond to the asymptotic behavior of a quark distribution 
amplitude oc (1 — a^) and for gluons oc (1 — a^)^, respectively. This correspondence is 
not demanded by theory. Therefore, we will consider both the cases, n = 1 and 2, for the 
construction of the gluon GPD. A parameterization of the t dependence of / is difficult. 
The multiplication of / as given in ()63p by a t-dependent form factor, although frequently 
used in default of a better idea, is unsatisfactory. Parameterizations of the GPDs ^71 IHH] 
as well as results from lattice QCD jHS] revealed that a factorization of / in /?, a on the 
one hand and in t on the other hand is most likely incorrect. Fortunately, the knowledge 
of the GPDs at t ~ suffices for our purposes as will become clear in Sect. IHl 

According to [32], the GPD is related to the double distribution by (since we will 
only work with GPDs at t ~ we omit the variable t in the GPDs in the following) 

H<^ix,0 = \<S>iO<x<0 r d(3 + Qi^<x<l) r d/3]^fi/3,a ""'^^ 



HD^{^)- (65) 
The definition of is completed by noting that it is an even function of x 

H^{-x,0 = H%x,i). (66) 
The integration limits in Eq. ()65|1 are given by 

x + i x-i x-i 

The limit Xi(x2) is the momentum fraction the emitted (reabsorbed) gluon carries with 
respect to the incoming (outgoing) proton. The last term in the definition is the 
so-called D-term jlU]. Its support is the region — ^ < a; < ^ and it ensures the correct 
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polynomiality property of the GPD. Since the D term is oc ^ and our interest hes in small 
skewness, we neglect it. 

We take the gluon distribution from the NLO CTEQ5M results jl^ and use an in- 
terpolation of it which has been proposed in Ref. jl2j and which is valid in the range 
Q2 = 4 GeV^ < Q2 < 40 GeV^ 

Pg{P) = r'^'i'^l-pfjZc,^''^ (68) 

where 

Co = 1.94, ci = -3.78 + 0.24 In (QVQo), ca = 6.79 - 2.13 In (QVQo) • (69) 

This is a very good approximation to the CTEQ gluon distribution for (3 < 0.5. At the 
largest value of we are going to use the interpolation (|UHj). namely 40 GeV^, and for 
f5 ~ 0.2 it deviates less than 5% from the CTEQ gluon distribution. For values of (3 in the 
range 10"^ - 10-\ the interpolation §^ agrees with the CTEQ gluon distribution within 
1%. In this region the Q^-dependence of 6 is approximatively given by 

5(Q2) = 0.17 + 0.07 In (Q^Ql) - 0.005 In^ (Q^Ql) . (70) 

The parameterization (jHHjl . (fHIIjl and (|7n|) effectively take into account the evolution of the 
gluon distribution in a large but finite range of as calculated in Ref. [Hj. At small (3 
the gluon distribution has a typical error of about 15% Within this error there is 

agreement with the analysis presented in Ref. |12I- An error assessment of the power 6 
provides an uncertainty of about 10 — 15% for it li^. 

For the various terms in the ansatz (j68|l the integrations occurring in Eq. (j65|) can be 
performed analytically, see Ref. jHEj- One finds 

'^^ 2^3 r(4 + V2-5) 

+ e(l - x)(2 + V2 - S) + } , X > e , 

+ (x^-x)}, x<e, (71) 

for the case n = 1. A similar but somewhat more complicated results is obtained for the 
case n = 2. This way we obtain an expansion of 

H3{x,0 = Ecn.i^™(x,0, (72) 

i 

with coefficients following from ()68j) . The evolution of the GPD is here approximated by 
that of the gluon distribution. The dominant contribution to vector meson electropro- 
duction comes from the imaginary part of the L ^ L amplitude (see Sect. ^ which is 
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Figure 3: Model results for the GPD in the small x range at t ~ and for the case 
n = 1. The sohd (dashed, dash-dotted) line represents the GPD at ^ = 5 (1 , 0.5) ■ 10~^ 
and at the scale 2 GeV. 



oc H^{^,^). Since, to a good approximation, H^{^,C,) equals a;Bj5'(a;Bj) at small ^ we have 
approximately taken into account the evolution. The full evolution of the gluonic GPD 
is complicated because of mixing with the flavor- singlet quark GPD. Its modelling would 
counteract any possible gain of accuracy obtained by the inclusion of the full evolution 
behavior. 

The GPD and its derivatives up to order n are continuous at x = ^. For ^ <^ x one 
can convince oneself that H^{x,C,) — > x g{x) up to corrections of order In the forward 
limit, ^ 0, the GPD reduces to the ordinary parton distribution xg{x). Results 
for in the case n = 1 are shown in Fig. El For x larger than ^ -C 1 there is practically 
no dependence on the skewness in contrast to the region x < ^ in accord with the general 
behavior of the model GPD just mentioned. The GPDs for n = 1 and 2 agree with each 
other on the percent level at small x. As we checked the numerical results for the cross 
sections obtained with both these GPDs are very similar; the differences in the imaginary 
(real) parts of the amplitudes are typically smaller than 1(7)%. In the following we will 
therefore show only numerical results for the case n = 1. 

Considering the coUinear limit of the subprocess amplitude (j45|) . one notices a double 
pole {x — C, + ie)""^ occurring in the T ^ T amplitude (j24j) [TB] . Partial integration leads 
to the integral 

dx d 



H^{x,Of{x,0 



(73) 



/o X — ^ + ie dx 

where / arises from the subprocess amplitude ()45|) . Since the derivatives of and / are 
continuous at x = ^ the integral exist. The transverse quark momenta are not needed for 
the regularization of the T — > T amplitude. 

A model for the GPD can be constructed analogously to (jHSl), the parton dis- 
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tribution g{P) is only to be replaced by its polarized counterpart Ag{P). The continuation 
to negative P is defined by 

Agi-(3) = AgiP) . (74) 
The GPD is antisymmetric in x 

H\-x,i) = -H\x,i). (75) 

We take A^f from Ref. [IH] and parameterize it analogously to (jHE)) 

(5Ag{P) = P'^'^'^il-P)'j2~c.p\ (76) 

where 

Co = 3.39 - 0.864 ln(gVg2) ^ ci = 1.73 + 0.2i IniQ^Ql) - 0.17 In^Q^Ql) 
C2 = 0.42-0.115 ln(QVQ2) _ 0.069 ln2(Q7Q2), (77) 

and 

5(q2) = 0.78 - 0.173 In(gVQg) . (78) 

The GPD if can than calculated analytically for either case, n = 1 and 2, with, for 
instance, the help of (f7T|). It is then represented by a sum analogously to (f7^. We 
finally remark that the polarized gluon distribution and hence is subject to much 
larger uncertainties than H^. 



6 Cross sections 

Vector meson electroproduction in the diffractive region has been extensively investigated 
at HERA tiaHnilllllllllSlini for large W and Q"^ but small xbj. Preliminary data from HI 
and ZEUS jl7||lHllini extend the range of for which electroproduction data are available. 
In order to confront the data with the theory developed in the preceding sections, one has 
either to extrapolate the data to t ~ or to take into account the t dependencies of the GPD 
and the subprocess amplitudes. The latter recipe is not straightforward. As we mentioned 
in Sect. EJit is not easy to find a plausible parameterization for the t dependence of the 
GPD because factorization in t and x, ^ most likely does not hold jHHllSni- We therefore use 
a variant of the first recipe and multiply the t ~ amplitudes (j37j) by exponentials 

~exp[t5f/2], (79) 

with slope parameters, BY {i = LL, LT, TT for L ^ L,T ^ L,T ^ T transitions, 
respectively), adjusted to experiment. The ansatz (f7^ is in accord with the expected 
exponential behavior of the GPDs ISZHSOl- Differences in the slope parameters arise from 
the t dependence of the subprocess amplitudes. 
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In the one-photon exchange approximation the ep epV cross section integrated over 
the azimuthal angle , reads 



d^a{ep epp) «eim I + {I - yf 



dydQ'^dt 



2tx 



dax dai 
'-e- 



dt 



dt 



where high-energy, small Xbj approximations have been applied to the phase space factor. 
Under the same kinematical conditions the ratio of longitudinal to transversal polarization 
of the virtual photon is given by 



i-y 



where y is the fraction of longitudinal electron momentum carried by the photon 



(81) 



y 



q-p 



(82) 



Here, kg is the momentum of the incident electron and s = {kg + pY . The 7*p Vp 
partial cross sections in (|Sm) for transversally and longitudinally polarized virtual photons 
are related to the amplitudes (j37j) by 



daT 
If 
dcTL 
dt 



1 



\M 



H 



(83) 



where we made use of Eqs. (j77|) . (PT|) and Terms of order {H^)^ have been neglected 
in the cross sections (j83p : there is no interference between the and contributions. 

The differential cross section data for ep epV exhibit a characteristic diffraction peak 
at small t. The slope of the diffraction peak is found to be nearly independent of W but is 
mildly varying with Q^. Most of the differential cross section data for p and production 
are compatible with a single exponential within errors. The combined HI and ZEUS data 
on the slopes in the range 4 GeV^ ^ ^ 40 GeV^ can be condensed into 



B 



LL 



7.5 GeV~2^1.2 GeV-^ln 



3.0 GeV' 



Q'^ + m 



V 



This parameterization is in rather good agreement with experiment, keeping in mind that 
the experimental slopes are not always extracted from cross section data in the same range 
of t. Possible deviations from a single exponential behavior of the cross sections then lead 
to different slopes. We naturally assign the slope (jH^ to the dominant L — >■ L transition 
amplitude. The slopes of the other amplitudes are not well determined as yet. A detailed 
analysis of the spin density matrix elements presented in Sect. seem to favor the choice 

These slope values lead to results from our GPD based 



BYl slightly. 
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Figure 4: The differential cross section ()80|) for ep epp versus —t integrated over tlie 
kinematical region available to the ZEUS experiment from which the data are taken [TH] . 
The solid line is our fit to the data at W = 75 GeV and = 6 GeV^ (see text). The 
dashed (dot-dashed, short-dashed) line represents the individual contributions from the 
L—>-L{T—>-T,T^L) amplitudes. 

approach in fair agreement with the HERA data. It is to be stressed that the magnitude 
of the transverse cross section is controlled by the product of parameters {fvT/MvY / B^j, 
leaving aside the mild dependence of the slope. The just described fit is based on the 
choice My = mv (see the remark subsequent to Eq. (jBl})- Taking a smaller value for My 
and a corresponding value for the decay constant, the slope B^j^ can be closer to that one 
for the L L amplitude. For instance, choosing My = mv/'2, one may use B^j, = B\^ 
( for fpT = 170 MeV and f^r = 190 MeV) and obtaines almost identical results for the 
cross sections. 

As a check of our choice of the slopes we show the ZEUS data ^3] for the differential 
cross section of p production in Fig. |3J These data indicate deviations from a single 
exponential behavior. They are integrated over the W and region accessible to ZEUS; 
W varies between 32 and 167 GeV in dependence on which varies between 3 and 
50 GeV^. The associated normalization uncertainty is of no bearing to us since we are 
interested in the process YP ~^ Vp. The forward amplitudes (j37|) evaluated from the 
model GPD shown in Fig. |2l multiplied with the exponentials (|79|) . lead to the results 
for the ep — > epp differential cross section shown in Fig. |3] The agreement between our 
result and experiment is not too good. Obviously, the value of the slope taken from Eq. (j5^ 
at = 6 GeV^, is a bit too small. However, the data shown in Fig. HI need confirmation. 
We can also see from the figure that our result although obtained with different slopes, do 
not deviate from a single exponential behavior substantially. Also shown in Fig. 0] are the 
three individual contributions L — L, T — T and T — > L separately. As expected the 
L ^ L contribution dominates. The T ^ T contributions amounts to about 25% of the 
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Figure 5: The integrated cross section for 7*p pp (left) and 7*p (pp (right) versus 
at ~ 75 GeV. Data taken from [T21 El (filled squares) and [TSl ES] (open squares) for 
p and production, respectively. The solid lines represent our results. 

L L one at t ^ 0. Due to the smaller slope B^j, takes the lead for —t larger than about 
0.4 GeV^. The T L contribution is shown only for comparison, it is of no account to 
the cross sections. 

Let us now turn to the discussion of the process 7*p Vp. The integrated cross section 
for this process is related to the integrated partial cross sections (jH^ by 



(t(7*p Vp) = aT{l*p Vp) + taL{l*P Vp) . (85) 

The HI Ha iU and ZEUS [131 iS] data on the cross sections for YP ^ pV {V = p, 0), 
integrated over the diffraction peak, are compared to our results in Fig. We repeat our 
results are evaluated from the handbag amplitude (|Ff|) multiplied by the exponentials (f?^ 
and using the GPD shown in Fig. |21 Good agreement between model and experiment 
is achieved for both processes provided is larger than about 4 GeV^. 

The HERA experiments also measured the decay angular distributions of the p and 
mesons and determined their spin density matrix elements. This information allows for a 
determination of the cross section ratio 

R(V) = (86) 

axi^y p — >■ Vp) 

from which, in combination with ()85p . the longitudinal cross section, ai, can be isolated 
as well. The HERA data for (Tl and R are compared to our results in Figs. IHl and [71 
respectively. Again reasonable agreement is to be observed for Q"^ larger than 4 GeV^. 
The ratio R increases with Q'^ since the transverse cross section is suppressed by as 
compared to the longitudinal one, see the hierarchy (j4H) . 

The experimental results on cross section ratio are derived from data on the spin density 
matrix element Tqq. The extracted ratio is therefore a ratio of the differential cross sections 
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Figure 6: The integrated cross section for longitudinally polarized photons versus Q2 at 
14^ ~ 75 GeV. Left: 'yip ^ pp; right: 72p — > 4>p. The data are taken from HI [T^liil liH] 
(filled squares) and ZEUS jlHl ES] (open squares), respectively. The solid hues represent 
our results. 




Figure 7: The ratio of longitudinal and transverse cross sections for p (left) and (right) 
production versus at ~ 75 GeV. Data taken from [121 EH EEl (filled squares) and 
[ISl ESI (open squares) for p and production, respectively. The open triangles represent 
preliminary ZEUS data |1H1 EH] for p and electroproduction. The solid (dashed) lines 
are our results for the ratio of differential (integrated) cross sections, R{R). The ratio R 
is evaluated at t = —0.15 GeV^. 



24 



which equals the ratio of integrated cross sections, i?, only if both the differential cross 
sections show the same t dependence ^. This is however not the case if the slopes differ. 
Therefore, R, measured at t ~ —0.15 GeV^, is about 10 — 20% larger than R. In Fig. Owe 
also display our prediction for R. Very good agreement with experiment is to be seen now. 
It is to be stressed that the uncertainties of the gluon distribution |1T] entail a typical error 
of about 30% for our predictions for the cross sections. In the ratios these errors cancel to 
some extent. As we remarked a fit with B^j, ~ B^i is also in agreement with the present 
data provided the value of the product {fvT/MyY/B'^j, is kept constant. The ratio R for 
this fit practically falls together with R in the fit presented above. 

In Fig. ISl we display an Argand diagram of the three forward amplitudes for p elec- 
troproduction at = 4 GeV^ t = -0.15 GeV^ and W = 75 GeV. Both, A^^,o+ and 
++ are dominantly imaginary while the T ^ L one is nearly real. The latter phase 
is a consequence of the branch point of a/x^ — in Eq. ()44p . The hierarchy (PT|) is here 
seen again. The phase of the p-production amplitude A4q^ at t ~ is shown in more 
detail on the right hand side of Fig. |H1 The real over imaginary part ratio increases with 

and takes values between 0.2 and 0.4 in the region of interest. The real part of the 
L ^ L amplitude therefore contributes only about 10% to the cross sections. 

A number of comments concerning the leading-twist contribution OE] are in order. As 
we mentioned above it is given by the collinear approximation of the dominant amplitude 
M-o^o^- The salient features of the leading-twist contribution are passed on to the full L 
L amplitude, the quark transverse momenta and Sudakov suppressions essentially affect 
its absolute value. The examination of the leading-twist contribution therefore elucidates 
many properties of our results in a simple way. Neglecting the k± terms in ()43p and using 
the standard definition of a meson distribution amplitude 

^""^ ^^vdr) = [^^vL{r,kl), (88) 



we obtain the subprocess amplitude q_,_ in collinear approximation from Eq. (ji^ and, 
inserting it into (HH), the leading- twist contribution to the L ^ L amplitude 



.coll ,1 /_\ ^ H3{x,C) 



The 1/r moment of the meson's distribution amplitude ^vl occurring now, is denoted by 
For the wavefunction (|T7jl the associated distribution amplitude is the asymptotic 
form PH|) which leads to a value of 3 for the 1/r moment. 

^For single exponentials the relation between R and R is given by 
R = Btt/Bll cxp [-{Bll - BTT)t] R{t). 
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Figure 8: Left: The p production amplitudes for the three transitions L ^ L, T ^ T and 
T ^ L at g2 = 4 GeV^ t = -0.15 GeV^ and W = 75 GeV. Right: Real over imaginary 
part of the amplitude A^o+,o+ ^^r p production versus at t ^ and W = 75 GeV. 

We can now easily understand the growth of the real over imaginary part ratio with 
Q^. Applying the derivative analyticity relation frequently but unjustifiedly termed 
the local dispersion relation [521, to the imaginary part of the leading-twist amplitude (|H^ 



M 



coll 

o+,o^ 



7T d 



-Xb] 



and using the low-^ behavior of the model GPD H^{C^,^) = co{2^)~^ (see Eq. fl7T|)). we find 

ReM^V/ ImMotW ^ ^vr ^(Q^) . (91) 

The increase of 6 with ( see Eq. (ffn|) ) which has been calculated by the CTEQ group 
[¥T] with the help of QCD evolution, is what we see at the right hand side of Fig. |H1 
Up to corrections from the real part the integrated longitudinal cross section reads 



= ^-^^^A^sfvLCv{l/r)y,]' \H%^,0\' , (92) 



in coUinear approximation. The ratio of the and p cross sections is given by {ftjiLC^IfpLCp)"^. 
Our results shown in Figs. |31 and El approach this value with increasing . Due to the 
behavior of H^{^,^) at small ^ the cross section behaves as 

af'ozW'^'^'K (93) 

at fixed and small xbj. The power behavior comes about as a consequence of the 
behavior of the GPD and the underlying gluon distribution. We note in passing that in the 
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Figure 9: Left:The integrated cross section for 7*p — > pp versus W for five values of Q^. 
Data taken from The solid lines represent our result. Right: The power 5 versus 
for p electroproduction. The solid curve represents the power as determined in Ref. jHj 
with an error estimate given by the shaded band. Data are taken from ^2] (•) and jT^ 



Regge picture [2] the exponent 5{Q'^) is associated with Pomeron exchange. The intercept 
of the Pomeron trajectory is related to 5 by ap(0) = 1 + 5{Q'^). In the Regge model 5 is 
a free parameter. 

In Fig. El we display the cross section for 7*p ^ as a function of W for sample 
values of . Fair agreement between experiment and our predictions is to be seen. The 
W dependence of the predictions from the full approach is very close to that given in 
(|93p . Deviations from the power law at lower values of W , to be observed in Fig. El arise 
from various corrections to the leading- twist contribution we take into account, such as 
the quark transverse momenta, the T — T amplitude and the real parts of the L ^ L 
amplitude. This interpretation of the power behavior of ctl is supported by a comparison 
of 5 as taken from the analysis presented in Ref. with the powers obtained from fits 
to the cross section data [121 CSl- Rough agreement between both the results is to be seen 
in Fig. El although the errors of the HERA data do not permit a definite conclusion as yet. 
Preliminary HERA data seem to improve the agreement. 

For very small ^ one can estimate the size of the collinear contribution using the leading 
terms in the model GPD ()71|). One obtains 

(94) 
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Figure 10: Ratios of the GPDs and the parton distribution for the n = 1 and 2 models at 
a scale of 4 GeV^. 

where cq = co/[(l — 5/3)(l — S/2)] is the coefficient of the first term in the power series of 
()72|) associated with Eq. (jUHj) . This cross section is rather large, well above experiment. 
The quark transverse momenta and the Sudakov factor suppress it such that agreement 
with experiment is found, see Fig. El 

Exploiting the leading ln(l/xBj) approximation, Brodsky et al. 0] found a result ^ that 
equals (jHH) except that H^{^,$,) is replaced by the usual gluon distribution XBjfi'(a;Bj) (see 
also f53^). At very small ^, i.e. if ^ is so small that the first terms in Eq. ()68|) and in 
the corresponding GPD (f7^ suffice; the usual gluon distribution and the GPD only differ 
by about 20% resulting from the difference between cq and Cq- For larger ^, however, the 
difference between both the functions becomes substantial, growing up to about a factor 
of 1.6 - 2 at ^ = 0.1, see Fig. HUl The use of the ln(l/xBj) approximation at values of 
^ around 0.1 may therefore lead to an underestimate of the gluonic contribution to cross 
sections by a factor 3 to 4. We repeat that, in contrast to the ln(l/xBj) approximation, 
we do not require small ^ in principal. We only restrict ourselves to small ^ in order to 
avoid complications with potential contributions from quarks emitted and reabsorbed by 
the proton. The enhancement effect apparent in Fig. ^1 is known as the skewing effect 
and has been discussed by several authors |Hl US IMl ESI- The size of the skewing effect 
estimated in these papers, is compatible with our model result for ^ ^ 1. 

In any case the leading- twist as well as the ln(l/xBj) result provide cross sections that 
are too large. In order to settle this problem for the ln(l/xBj) approximation Frankfurt et al 
^ estimated a correction factor by allowing for quark transverse momenta in the loop. This 
mechanism bears resemblance to our approach as we mentioned in Sect. El The correction 

^Note, in 0] the decay constant includes the flavor weight factor Cy- 
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factor obtained in [H] is large enough to achieve agreement with experiment. This factor 
has also been used by Mankiewicz et al. ^3] in an explorative study of (Tl in an otherwise 
collinear GPD approach. Martin et al. J8j exploited the ln(l/xBj) approximation in their 
analysis of vector-meson electroproduction by including parton transverse momenta and 
an unintegrated gluon distribution. 



7 Spin density matrix elements 

With the help of Eqs. (j2ZI), PT|) and the spin density matrix elements extracted 
from the decay angular distributions measured with unpolarized leptons and protons 
simplify to {R is defined in flH7|) ) 
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while 



.04 
1-1 



11 



11 



rl_^ = Im rl_-^ = 



(96) 

because of the neglect of L — > T and T —>■ —T transitions. The relations ()95|) . obtained in 
the GPD approach under the assumption of the dominance of the terms, coincide with 
those found assuming dominance of natural parity t-channel exchanges and the neglect 
of proton helicity flip ^21 IHl] • The contributions from enter the spin density matrix 
elements only as bilinears, there are no interferences between and terms. 

The data for the spin density matrix elements from HI ^21 ISj and ZEUS ^1] , are 
shown in Figs. ^2 and ^1 and compared to the results from the GPD based approach. 
The general pattern of the data is reproduced. The dominance of the L ^ L transition 
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amplitude is clearly visible in the angular distribution of the production and decay of the 
vector mesons, in particular in the value of Tqq which tends towards 1 with increasing . 
This behavior is well reproduced by our approach as we already discussed in connection 
with the cross section ratio R. 

The T L amplitude is probed by the matrix elements Tqq and Tqq. While the first 
matrix element is approximately oc ++P/|A^o+ o+P' ratio Im ++/Im q_i_ 

essentially controls the second. Both Tqq and ItqqI are found to be rather small. The ratio 
of the T T and L L amplitudes is approximately measured by rl_^ and Re r^Q, 
quadratically in the first case, linearly is the second one since both the amplitudes have 
about the same phase as is shown in Fig. |H1 The fair agreement between theory and 
experiment for these spin density matrix elements tells us that our approach provides the 
correct sizes and relative phases of the T ^ T and L L amplitudes. The matrix 
elements Re r^g = —Re r^Q = Im rf^ measure an interference term between the T ^ T and 
T ^ L amplitudes which is very small. Also this prediction is in acceptable agreement 
with experiment. 

The t dependence of the spin density matrix elements confirm the above observations, 
see Fig. ^| The T —>■ L sensitive matrix elements behave proportional to or t while 
those controlled by ratios of the T T and L L amplitudes exhibit an t dependence 
according to the different slopes chosen for them. As we mentioned in Sect. IHlthe freedom 
in choosing a suitable value of My also allows fits with B^j, ~ BYl- While the transverse 
cross section is nearly insensitive to this choice provided the product {fvT/MvY / Bj.rp 
is approximately kept fixed, does the t dependence of some of the spin density matrix 
elements (e.g. Tqq, rl_^) change; they become very flat in t. Given the accuracy of the 
present data |.12j such a behaviour is not in conflict with experiment. 

Finally, in Fig. we show our predictions for electroproduction at = 10 GeV 
characteristic of the COMPASS experiment. 

Other theoretical analyses [71 EHl EZl HH] of the spin density matrix elements base on 
variants of the ln(l/a;Bj) approximation. The variants differ from each other in the detailed 
treatment of the subprocess '-f*g —>■ qqg. The same hierarchy of the amplitudes are obtained 
as we find and, in general, rather similar results are obtained for the spin density matrix 
elements. Worth mentioning is the different phase of the T —>■ L amplitude and a somewhat 
different t dependence of the matrix elements. 



8 The helicity correlation 

Last we want discuss the role of the GPD H^. For this purpose we consider the initial 
state helicity correlation An which can be measured with longitudinally polarized beam 
and target. After integration over the azimuthal angle this correlation reads 



A^^[ep ^ epV] ^^^^ 1-^++,++! + l-Mo+,++r - - l-Mo+,-+ 



(97) 
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Figure 11: The spin density matrix elements of electroproduced p mesons versus 

at 

W ~ 75 GeV and t ~ -0.15 GeV^. Data, taken from [12] (filled circles) and 14^ (open 
circles), are compared to our results (solid line). Preliminary data on rgg from ZEUS jlH] 
(open triangles) are also shown. 
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Figure 12: The spin density matrix elements for electroproduction versus at W ^ 
75 GeV and t ^ -0.15 GeVl The HI data il are compared to our results (solid line). 
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Figure 13: The spin density matrix elements of electroproduced p mesons versus t at 
Q2 = 5 Qgv2 and W 75 GeV. Data taken from Hg (filled circles). The solid (dashed) 
lines represent our results for the choice B^j, = [BYj^)- 
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Figure 14: The spin density matrix elements of electroproduced mesons versus at 
~ 10 GeV, y ~ 0.6 and t ~ —0.15 GeV^. The sohd hues represent our results. 
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Figure 15: Left: The helicity correlation An for p electroproduction versus at W = 
15 GeV, t ~ and y ~ 0.6. Data taken from SMC [52]. Right: An for production at 
W = 5 GeV (sohd hne) and W = 10 GeV (dashed hne); y ~ 0.6. The shaded bands reflect 
the uncertainties in our predictions due to the error in the polarized gluon distribution 



where the amplitudes and cross sections refer to the process 7*p — > Vp and are given in 
Eqs. and (jHHI) . As can easily be seen from Eq. An = if the terms are 
neglected as we did in the preceding sections. Yet in contrast to the cross sections and 
spin density matrix elements where the correction are bilinear in the terms and, hence, 
extremely small, the leading term in A^l is an interference between the and the 
terms. In fact, with the help of Eqs. ()27|) and ()58|) . one obtains from Eq. ()97p 



Re 



All [ep ^ epV] = 2^1^ , i ^ i . .IT . • (9' 



Obviously, this ratio is of order {kj_)/Q'^ (H^) / (H^) and, therefore, very small values for 
All are to be expected. Indeed exploiting the model GPDs presented in Sect.Elwe conflrm 
this assertion as can be seen from Fig. fTHl where results for All for p and electroproduction 
at t ~ are displayed. The results for p production, only shown at W = 15 GeV, is 
compared to the SMC data jSlI. At this energy and in the range of shown in the plot, 
the contribution from the quark GPD is expected to be small ^Tj. Our results for All 
are not in disagreement with experiment given the admittedly large experimental errors 
and the rather large value of the skewness. Results for electroproduction are shown 
at energies typical for the HERMES and COMPASS experiments. The dominance of the 
gluon over the sea quarks permits this. At IV = 5 GeV All is not very small since the 
major contribution to it comes from the region 0.1 where Ag/g is not small. 

The proton helicity flip contribution, related to the GPD E^, may change these results 
but likely not substantially. The helicity correlation will increase with growing momentum 
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transfer if the slope of the T ^ T amphtude is smaller than that of the L ^ L one. Besides 
allowing predictions for All this calculation also supports our assumption of negligible 
contributions from to cross sections and spin density matrix elements. 

9 Summary 

We analyzed electroproduction of light vector mesons at small xbj within a GPD based 
approach. In this kinematical domain the gluonic GPD H^, parameterizing the response 
of the proton to the emission and reabsorption of gluons, controls the process. The gluonic 
GPD, not calculable at present, is constructed from an ansatz for the double distributions 
currently is use. In order to examine the influence of the model GPD on the numerical 
results for vector meson electroproduction we used two different versions for it (n = 1 and 
2). The differences in the numerical results obtained from these two models, are on the 
percent level. The subprocess amplitudes for Yd ~^ ^9 ^-^e calculated by us to lowest 
order of perturbative QCD but transverse momenta of the quark and antiquark that form 
the vector meson, are taken into account as well as Sudakov suppression which sum up 
gluonic radiative corrections. 

The GPD approach reproduce all main features of vector meson electroproduction at 
small XBj known from phenomenology. The dominance of the contributions from the GPD 

over those from and , leads to the relations ()58|) and, hence, to results equiva- 
lent to those obtained assuming the dominance of natural parity exchange. Approximate 
s-channel helicity conservation holds due to the hierarchy ()4H) the amplitudes respect in 
our GPD based approach. The behavior of the longitudinal cross section as a depen- 
dent power of W at fixed is a consequence of low ^ properties of the GPD and QCD 
evolution. The numerical results we obtain from our approach are in reasonable agreement 
with the small xbj data on cross sections and spin density matrix elements for electropro- 
duction of p and (f) mesons measured by HI and ZEUS. The t dependence of vector meson 
electroproduction is not yet satisfactorily settled. In principle it is generated by a combi- 
nation of the t dependence of the GPD and, with lesser importance, that of the subprocess 
amplitudes. Due to the lack of a plausible parameterization of the t dependence of the 
GPD we have evaluated the electroproduction amplitudes at t ~ and multiplied them by 
exponentials in t. Improvements on this recipe are demanded and will be unavoidable as 
soon as detailed differential cross section data are at hand. 

We also compared in some detail our approach to the leading-twist contribution and 
to the leading ln(l/xBj) approximation. The latter is rather close to the GPD approach 
at low xbj and small t but not identical. For xbj larger than about 0.01 the replacement 
of by 2C,g{2^) becomes inappropriate. The GPD has, in contrast to the leading 

ln(l/a;Bj) approximation, the potential to investigate the t dependence of electroproduction. 
The lack of understanding of the GPD's t dependence prevents this at present. 



36 



Acknowledgements 



We thank A. Borissov, A. Efremov, O. Nachtmann, H. Meyer, E. Paul, A. Sandacz and L. 
Szymanowski for discussions and A. Bruni for providing us with prehminary ZEUS data. 
This work has been supported in part by the Russian Foundation for Basic Research, Grant 
03-02-16816, the Integrated Infrastructure Initiative "Hadron Physics" of the European 
Union, contract No. 506078 and by the Heisenberg-Landau program. 



References 

D. Schildknecht, G.A. Schuler and B. Surrow, Phys. Lett. B 449, 328 (1999). 

A. Donnachie and P. V. Landshoff, Phys. Lett. B 437, 408 (1998); J. Bartels and H. 
Kowalski, Eur. Phys. J. C 19, 693 (2001), |h^h/0010345| 



A. Donnachie and P. V. Landshoff, Phys. Lett. B 185, 403 (1987); O. Nachtmann, 
contribution to the Ringberg workshop on HERA physics 2003, |hep-ph/0312279[ 

S. J. Brodsky, L. Frankfurt, J. F. Gunion, A. H. Mueller and M. Strikman, Phys. Rev. 
D 50, 3134 (1994) |hep-ph / 9402283]. 

B. Z. Kopeliovich et al, Phys. Lett. B 309, 179 (1993); J. Nemchik et al, Z. Phys. C 
75,71 (1997); I. Royen and J.R. Cudell, Nucl. Phys. B 545, 505 (1999). 

L. Frankfurt, W. Koepf and M. Strikman, Phys. Rev. D 54, 3194 (1996) 
|hep-ph/9509311j. 



D. Y. Ivanov and R. Kirschner, Phys. Rev. D 58, 114026 (1998) |h ep-ph/9807324 



A. D. Martin, M. G. Ryskin and T. Teubner, Phys. Rev. D 62 (2000) 014022 
[hep-ph/9912551j. 

A.V. Radyush kin, Phys. Lett. B 385, 333 (1996) |hep-ph/960543l| ; Phys. Rev. D 56, 
5524 (1997) |hep-p h/9704207|. 

D. Miiller, D. Robaschik, B. Geyer, F. M. Dittes and J. Horejsi, Fortsch. Phys. 42, 
101 (1994) |l|hep-ph/9812448|; X. Ji, Phys. Rev. D 55 7114 (1997) ^-ph/960938l|. 

J.C. Colli ns, L. Frankfurt and M. Strikman, Phys. Rev. D 56, 2982 (1997) 
|hep-ph/96ll433). 



C. Adloff et al. [HI Collaboration], Eur. Phys. J. CIS, 371 (2000) | hep-ex/9902019j . 

J. Breitweg et al. [ZEUS Collaboration], Eur. Phys. J. C 6, 603 (1999) 
|hep-ex/9808020| . 



37 



[14] Breitweg et al. [ZEUS Collaborations], Eur. Phys. J. C 12, 393 (2000) 
|hep-ex/9908026|. 

[15 
[16 



[17; 

[18 
[19 
[20 

[21 
[22 
[23 

[24 

[25; 

[26; 

[27; 
[28 

[29 

[30 

[31 

[32 



L. Mankiewicz, G. Filler and T. Weigl, Eur. Phys. J. C 5, 119 (1998) |' hep-pti/97H"227| . 

L. Mankiewicz and G. Filler, Phys. Rev. D 61, 074013 (2000) |hep-ph/9905287|. 
1. V. Anikin and O. V. Teryaev, Phys. Lett. B 554, 51 (2003) | hep-ph/0211028. . 

M. Diehl and A. V. Vinnikov, |hep-ph/0412"T62| 

H. W. Huang and P. Kroll, Eur. Phys. J. C 17, 423 (2000), |hep-ph/ 0005318 . 
J.B. Kogut and D.E. Soper, Phys. Rev. D 1, 2901 (1970). 

M. Diehl, T. Feldmann, R. Jakob and P. Kroll, Nucl. Phys. B 596, 33 (2001), 
[Erratum-ibid. B 605, 647 (2001)] ,hep-ph/0009255| . 



M. Diehl, Eur. Phys. J. C 19, 485 (2001) ]hep-ph/01 01335| . 



M. Diehl, Phys. Rept. 388, 41 (2003) |hep-ph/0307382.,. 



J. Botts and G. Sterman, Nucl. Phys. B 325, 62 (1989); H. Li and G. Sterman, Nucl. 
Phys. B 381, 129 (1992). 



R. Jakob and P. Kroll, Phys. Lett. B 315, 463 (1993) |hep-ph /9306259l; Erratum-ibid. 
B 319, 545 (1993). 

P.A.M. Dirac, Rev. Mod. Phys. 21, 392 (1949); H. Leutwyler and J. Stern, Ann. Phys. 
(N.Y.) 112, 94 (1978). 

V. L. Chernyak and A. R. Zhitnitsky, Phys. Rept. 112, 173 (1984). P. Ball and 
V. M. Braun, Phys. Rev. D 54, 2182 (1996) |h'ep-ph /9602323| . 

J. Bolz, J. Korner and P. Kroll, Z. Phys. A 350, 145 (1994). 

B. Guberina, J.H. Kiihn, R.D. Peccei and R. Riickl, Nucl. Phys. B 174, 317 (1980); 
S. Balk, J.G. Korner, G. Thompson and F. Hussain, Z. Phys. C 59, 283 (1993). 

S.J. Brodsky and G.P. Lepage, in: Perturbative Quantum Chromodynamics, ed. A. 
Mueller (World Scientific, Singapore 1989). 

M. Vanderhaeghen, P. A. Guichon and M. Guidal, Phys. Rev. D60, 094017 (1999) 
[hep-ph/ 9905372|. 

J. C. Collins and D. E. Soper, Nucl. Phys. B 193, 381 (1981) [Erratum-ibid. B 213, 
545; "Parton Distribution And Decay Functions," Nucl. Phys. B 194, 445 (1982). 

P. Kroll and M. Raulfs, Phys. Lett. B 387, 848 (1996) [hep-ph/9605264| . 



38 



[33] M. Neubert and B. Stech, in Heavy Flavours II, edited by A.J. Buras and M. Lindner 
(World Scientific, Singapore, 1997). 

[34] K. Schilling and G. Wolf, Nucl. Phys. B 61, 381 (1973). 

[35] M. Dahm, R. Jakob and P. KroU, Z. Phys. C 68, 595 (1995) |hep-ph/9503"4T8| . 

[36] I. V. M usatov and A. V. Radyushkin, Phys. Rev. D 61, 074027 (2000) 
|hep-ph/ 9905376|. 

[37] M. Diehl, T. Fe ldmann, R. Jakob and P. KroU, Eur. Phys. J. C 39, 1 (2005) 
|hep-ph/0408173 . 

[38] M. Diehl, T. Fe ldmann, R. Jakob and P. KroU, Eur. Phys. J. C 8, 409 (1999) 
||hep-p h/9811253 . 

[39] J. W. Negele et al, Nucl. Phys. Proc. Suppl. 128, 170 (2004) |h ep-lat/0404005j . 

[40] M. V. Polyakov and C. Weiss, Phys. Rev. D 60, 114017 (1999) !hep-ph/9902"45l| . 

[41] J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky and W. K. Tung, JHEP 
0207, 012 (2002) |hep-ph/0201195]. 

[42] A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Tho rne, Phys. Lett. B 531, 216 
(2002) hep-ph/0201127|; Eur. Phys. J. C 28, 455 (2003) |hep-ph/0211080l . 

[43] J. Bliimlein and H. B5ttcher, Nucl. Phys. B 636, 225 (2002) |hep-ph/0203T55| . 

[44] C. Adloff et al. [HI Collaboration], Phys. Lett. B 483, 360 (2000) |hep- ex/00050To|. 

[45] M. Derrick et al. [ZEUS Collaboration], Phys. Lett. B 380, 220 (1996) 
|hep-ex/9604008 . 

[46] S. Aid et al. [HI Collaboration], Nucl. Phys. B 468, 3 (1996) hep-ex/9602007 . 

[47] HI collaboration, contribution to the International Europhysics Conference on High 
Energy Physics, EPS03, Aachen (Germany) 2003. 

[48] ZEUS collaboration, contribution to the International Europhysics Conference on High 
Energy Physics, Budapest (Hungary) 2001 (Abstract 594). 

[49] ZEUS collaboration, contribution to the International Conference on High Energy 
Physics, Beijing (China) 2004 (Abstract 6-0248). 

[50] K. Goeke, M. V. Polyakov and M. Vanderhaeghen, Prog. Part. Nucl. Phys. 47, 401 
(2001) hep-ph/0106012|. 

[51] J. B. Bronzan, G. L. Kane and U. P. Sukhatme, Phys. Lett. B 49, 272 (1974). 



39 



[52] G. K. Eichmann and J. Dronkers, Phys. Lett. B 52, 428 (1974). 
[53] J. Bartels and M. Loewe, Z. Phys. C 12, 263 (1982). 

[54] A. D. Martin and M. G. Ryskin, Phys. Rev. D 57, 6692 (1998) [hep-ph/9711371]. 
A. G. Shuvaev, K. J. Golec-Biernat, A. D. Martin and M. G. Ryskin, Phys. Rev. D 
60, 014015 (1999) | hep-ph/9902410. . 

[55] A. V. Behtsky, D. Miiller and A. Kirchner, Nucl. Phys. B 629, 323 (2002) 
^hep-ph/0112108j . 

[56] J. Nemchik, N. N. Nikolaev and B. G. Zakharov, Phys. Lett. B 341, 228 (1994) 
[hep-ph/9405355|. 

[57] L Royen, Phys. Lett. B 513, 337 (2001) |hep-ph/0006044| . 

[58] L Ivanov, PhD thesis, hep-ph/0303053[ 

[59] A. Tripet [Spin muon collaboration], Nucl. Phys. Proc. Suppl. 79, 529 (1999) 
[hep-ex/9906008 . 



40 



